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Abstract. We propose an alternative perturbative expansion for QCD. All 
scheme and scale dependence is reduced to one free parameter. Fixing this 
parameter with a fastest apparent convergence criterion gives sensible results 
in the whole energy region. We apply the expansion to the calculation of the 
zero flavor triple gluon vertex, the quark gluon vertex, the gluon propagator 
and the ghost propagator. A qualitative agreement with the corresponding 
lattice results is found. 
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1 Introduction 



Perturbation theory is by far the most successful tool to get quantitative 
predictions from a field theory. Unfortunately, the results depend on the 
renormalization scale and scheme and the number of free parameters de- 
scribing this dependence grows with the order of truncation. In most cases 
one does not bother too much about this dependence and simply chooses 
a scheme (MS, MOM,. . .) which is supposed to give good results. More 
sophisticated approaches select a different ideal scheme for each perturba- 
tive series. Most cited in this context are Stevenson's principle of minimal 
sensitivity |jj and Grunberg's method of effective charges Q. 

In this paper we will reorganize the conventional perturbation series in 
a new alternative expansion, which has only one redundant free parameter, 
the expansion parameter itself. This parameter will be fixed by a fastest 
apparent convergence criterion (face). We find similar results as for the 
ordinary QCD perturbation theory in the UV-region and for intermediate 
energies. In the IR-region on the contrary, our expansion still gives sensi- 
ble results whereas the normal perturbation theory becomes useless if one 
approaches the Landau-pole. 

In section § we will first review some of the major aspects of the ordi- 
nary perturbation theory and then use this as starting point for the alter- 
native expansion. Some results will be presented in section ||. We perform 
calculations on the triple gluon vertex, the quark gluon vertex, the gluon 
propagator and the ghost propagator. Our results are compared with the 
conventional perturbation theory and with the corresponding lattice results. 
We also show that the reason for the IR finite results lies in a peculiar be- 
havior of the running expansion parameter y, resulting from the face. This 
behavior should be contrasted with the running coupling a that Shirkov 
and Solovtsov obtain by imposing analyticity Q. While we find a universal 
power behavior for y, they find an IR-finite value for a at zero momentum. 

2 Rewriting perturbation theory 

In the following we will consider the perturbative calculation of a renormal- 
ization scheme and scale invariant quantity 1Z, that is function of only one 
external scale q 2 , in a massless version of QCD . In ordinary perturbation 
theory one finds a row of approximations lZ n for 1Z, with 

n n = h (n) N ^ + ri ( q 2) h (n) + ^(^(n) 2 + _ _ _ + rn ( g 2)ft(») n ) . (2.1) 
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With N depending on the calculated quantity. The coupling constant 
is the solution of 

/3olng = £ + f MM) + 

\ []_ _Ai Po \ (99 , 

o W Pox Pox^ + p lX ^ + ... + p n x^)' [ -> 

with fi 2 -^h = — (Poh 2 + Pih 3 + P2h 4 + ...). As mentioned in the in- 
troduction, every truncation lZ n is highly scale and scheme dependent. 
One can for instance describe this dependence with the free parameters 
0] 0o In @2i /?3j ■ • • j Pn- The dependence of the coefficients on each of 
these parameters, cancels the dependence of the coupling constant up 
to order /j( n ) n+1 . 

In many cases one simply chooses a scheme (MS, MOM, . . .) and sets 
the scale /j, 2 equal to the external scale q 2 . To get reliable results, one must 
hope that the first, second, or third order approximation lies close enough 
to the exact result. The working hypothesis of perturbation theory is that 
the perturbation series is asymptotic to the exact result 1Z Q. One then 
obtains an error estimation, by assuming that 

= A n . (2.3) 

Other approaches fix the scale/scheme by imposing a condition on the 
truncated series. The minimal sensitivity condition [|J gives a different 
scale/scheme for every approximation. The method of the effective charges 
Hi, sets q 2 = fi 2 , the free parameters are then obtained by demanding the 
coefficients rj to be zero. One now also finds A n to give a good estimation 
of the error || . 

As a consequence of asymptotic freedom, conventional perturbation the- 
ory works well in the UV-region, which is reflected in low values of A n 
(n=l,2,3) for high values of q 2 . Unfortunately A n gets larger if we lower 
the external scale q 2 . This signals that one has to add non-perturbative 
power corrections to the conventional perturbation theory For inter- 
mediate energies, the sum rules || successfully relate many of these power 
corrections to a few condensates. A further lowering of q 2 towards the IR- 
region is catastrophic in most cases 0, A n diverges together with TZ n as one 
encounters the Landau-pole and perturbation theory becomes useless. 

2 Exceptions can be found in ||, where the minimal sensitivity criterion selects a scheme 
with an IR fixed point. 
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The alternative expansion we propose in this paper has no Landau-pole 
problem and gives sensible results up to q 2 = with a reasonable error 
estimation. We simply have to exchange h for a new expansion parameter 
y, defined by 



^0 ln To = ~ + IT ln (A)2/)- 
V Po 



From (|2.1|),( p^2|) and ( |2.4j ) we find after some calculations 



(2.4) 



K 



y N i + y 



Ai + Nk 



(h 
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+k((N + 1)^ + N^) +k 2 ^(N + i: 



,N 



+y 3 A 3 + ^i(iV + l)( 



1\2 



+k(A 1 (N + l)|l + (iV + 2)A 2 + iV(iV + 2)^ 
v Po Po 



--(-r) + - , 



-JV(iV + l)A 2 ) 



+A; 2 (^(iV + 2)(iV + l)+iV(Ar + ^| 
+k 3 [^(N + 2)(N + l)) +... 



(2.5) 



with 



% 2 ,2/) 



i + ^ln(/3 y)-/5 ln-^ 

y A A- 



MS 



(2.6) 
(2.7) 



The /^-coefficients (32,03, ■ ■ • are also in the MS'-scheme. 

We now find a row of approximations TZ n (y) with all the redundant 

2 

dependence residing in one single free parameter (3q ln ^ or equivalently y. 
All the other scheme dependence has disappeared, because it was eliminated 
from the expansion parameter. This might seem strange, since we explicitly 
refer to the MS scheme, but one can show that another choice for the 
reference scheme changes the coefficients and /3j in such a way that it 
exactly compensates the shift of k (k^k' = k + 2(3$ ln , A ). 



The scheme/scale dependence of ri_,T2,rz was derived in 



4 



To fix the expansion parameter y for each value of q 2 , we will have to 
impose some condition. We choose a face: for each approximation lZ n (y) we 
will set y equal to y n , the expansion parameter that minimizes the relative 
correction to the first order truncation of the series: 



mm 



K n (y)-K l {y) 



n n {y n )-n\y n ) 



(2- 



We will use the same error estimation as for the conventional perturbation 
theory: 



A n = 



U n+ \y 



n+l) 



n n {y ri 



Tl n + l {y 



(2.9) 



(Notice the difference: the fixing of y is done on the series, while the error is 
estimated on the results obtained after fixing.) As for any other possible fix- 
ing condition (e.g. minimal sensitivity, another face,...), there is no rigorous 
mathematical motivation for the condition we use. The true motivation lies 
in the fact that it generates sensible results with a good error estimation. 
This is the case for all the calculations we have performed so far. If we 
use, for example, the more obvious face, where one minimizes the relative 
correction to the zero order truncation y N , we do not find sensible results 
for the whole range of energies. There is a discontinuity at the point that 
separates the high energy region where this correction minimizes to zero, 
and the low energy region where it minimizes to a nonzero value. Let us 
now present some results obtained from (|2.5| ) and ([2. §|). 



3 Some results 

We now demonstrate the y-expansion on some quantities that also have been 
calculated on the lattice. All the needed two and three loop results have 
been calculated by Chetyrkin and Retey || . Everything is in Landau gauge 
for N c = 3 and Nf = 0. We will take the method of effective charges to 
be exemplary for the ordinary perturbation theory, but similar results are 
found with any other approach to the conventional perturbation theory. 



3.1 Triple gluon vertex 

There are several ways in which one can associate a (renormalization) scale 
and scheme invariant coupling constant with the triple gluon 3-point func- 
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tion 

G$? C ( P , q) = i 2 jdxdye-^+w) < T[A^(x)Al(y)A c p (0)] >, (3.10) 

or more precisely with its related vertex function (p, q, —p — q), defined 
by 

G$tM = D%,(-p)D b v l,(-q)D c p f p ,(-p - q)T% p ,{p,q, -p - q), (3.11) 
where 

D*(q) = i Jdxe^ < T[A^(x)A b v (0)] > . (3.12) 

If one sets one external momentum to zero, one finds |§ that the vertex- 
function can be written as: 

Yf° p {q, -q, 0) = -igf abc [(2g^q p - 9fip q v - g^T^q 2 )) 

-(9*> ~ ^rhpUq 2 )) (3.13) 

The coupling that was calculated on the lattice || [[(J is found to be section 
6.4): 

a s tf) ee ^OMgg^ 2 ) = h^-q 2 ) - ±T 2 (-<?)) 2 Z(-q 2 ) s , (3.14) 
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where h = i g _ 2 , (3.15) 



D*(q) = 5 ab (9»v- q -^)^p. (3.16) 



16vr 

[Uuv 

q q 

One can easily check the scheme and scale independence of a s . The three- 
loop result for /jMOMgg - n ^ e MS-scheme for /x 2 = q 2 is,||: 

fc MOMgg = h + h2[ 70 s^^^ 

6 bib 4 

l4r 304676635 299961 81825 , , , mS 

M ^9T2 64- C3 -^r C5] ' (3 ' 17) 

where Q is the Riemann zeta function. From this we can read of the coef- 



ficients A\ , A2 , As that are needed in ( [2.5D . The /3-coefficients in the MS- 
scheme have been calculated up to four loops in [11]: 

2857 n 149753 . 
O = 11, fa = 102, fh = — , 03 = — p- + 3564C 3 . 3.18 

2 6 
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We will compare our two and three loop results for a s {q 2 ) obtained from 
( |2.5[) (with N = 1) and (|2,8| ) with the results obtained from the method of 
effective charges |2| or equivalently in the MOMgg-scheme defined on the 
triple gluon vertex. The two and three loop MOM-scheme results are found 
as the solution of ( |2.2[ ) with n = 2, 3. The MOM-scheme /^-coefficients can 
be easily obtained from ( |3.17 ) and ( 3.18 ): 



.MOMgg 186747 1683 
lh ~ 6 4 4 



.MOMgg 20783939 1300563^ 900075 . 
* = 32— C3 32~ C5 - (3 - 19) 



While the A parameter is given by |j~2] ]: 

2p ln\™^ = ™ (3.20) 

Our two and three loop results are plotted together with the two and three 
loop MOiM-results in figs. | and g. 

We can clearly distinct three regions. For q > 30A^jg one finds the 
UV-region: the four results for a s coincide and the perturbation theory is 
completely reliable. The intermediate energies region goes from q ~ 30A^jg 
down to q ~ lOA^jg. A difference grows between the two and three loop 
results, but for both orders the y-expansion results still coincide with the 
MOM results. Power corrections are expected. For q < lOA^jg we find 

ourselves in the IR-region. The MOM results diverge while the y-expansion 
results continue to behave in a sensible way. 

The same conclusions can be read from fig. |3], where A 2 (see ( |2.3| )) is 
plotted, both for the MOM scheme and for the y-expansion. In the IR- 
region the error estimation diverges for the MOM scheme, while it stays in 
an acceptable interval for the y-expansion. 



It is the face (2^) that keeps the error estimation under control in the 
IR. This criterion, and in fact every other sensible criterion, will select for 
each momentum q a value for y that makes the higher order (n > 1) terms 
in the series (^1)|) as small as possible. Both for small values of q (q <C Aj^Jg) 
and for large values (q 3> Aj^g) it is the value of 

yk{q\ y) = l + y ^ ln(/? y) - y(3 In , (3.21) 

that determines the size of these higher order terms. For q S> ^j^g the large 
logarithm will be compensated by the y that multiplies it. One finds the 
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Figure 1: a s (q) (q in units of Ajjjg)> for 2 and 3 loops in the y-expansion 
and in the MOM-scheme. 
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q(A) 

Figure 2: Zooming in on the intermediate energy region of fig. [l] 
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Figure 3: q (in units of A^jg) — > A 
y-expansion and the MOM scheme. 



„(2) „,(3) 



7v 



(in %), for both the 



usual high energy running of the expansion parameter: 

1 



y = — 3 , (3-22) 



with c a constant. This gives 



% 2 , v)y q «°° — V- ( c " I ln < ln F") ) • ( 3 - 23 ) 

MS 

For g -C A ^jg the same cancellation can not occur since y must be positive, 
the large logarithm will now be compensated by the logarithm in the second 
term, we find a power behavior for y: 

V^A-izf (3-24) 

\ MS r 

and 

k(q 2 ,y)y q ^° 1. (3.25) 
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Figure 4: y(q) (q in units of Aj^g) for the tree loop truncation. Also de- 
picted: the low and high energy fits ( |3.24| ) and ( p. 22 ). 



This high and low energy behavior of the expansion parameter y is com- 
pletely universal, it is independent of the order of truncation, of the coeffi- 
cients Ai and to a certain extent of the criterion that was used. The running 
of y is depicted in fig. |] together with the fitted low ( |3.24|) and high ( |3.22|) 
energy behavior for the three loop truncation. 

If one would use the series itself to estimate the truncation error, ( |3.25|) 
would seem to invalidate the expansion for low energies, since the higher 
order terms become order 1. However, if you look at the the row of trun- 
cations ( |2.9D to estimate the error , the expansion remains valid (at least 
for low orders) since 2.5% < A 2 < 7.5% (see fig. |2.3| ). We have found a 
same behavior of A 2 for every other possible vertex-coupling that could be 
calculated from ||. 

We will finally compare our results with the lattice results of Q. This 
requires a fit of A-^y^, which was done for the two and three loop MOM 
results in || and || in the intermediate energy region (3 Gev-10 Gev). It 
was found that the MOM-results could be best fitted to the lattice results 
if a power correction was added. The fitted two and three loop values of 
A-tj-j are: 235 Mev and 238 Mev. The three loop power correction is 30% 
less then the two loop one. 

Since in the intermediate energy region the results of the y-expansion are 
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the same as the MOM results we can rely on the aforementioned fits. We 
will use the same value A-^g^=237 Mev, for every order. The two and three 
loop results of the y-expansion are plotted together with the lattice results 
in fig. [|. As expected, the difference between our results and the lattice 
result can be fitted as a power correction for q > 3 Gev. The amplitude 
of our maximum is significantly smaller then the amplitude for the lattice 
maximum in the IR region. But both maxima seem to approach each other, 
our amplitude grows larger with the order of truncation while the lattice 
amplitude becomes smaller for larger volumes (smaller P). 

3.2 The quark gluon vertex 

Again, there are several ways one can associate a scale and scheme invariant 
running coupling with the (zero flavor^) quark-gluon vertex A®y, which is 
defined by: 

G^jp, q) = Swi-P^ 'fip, q, ~q - p)S J , J (q)D a ^(p + q) , (3.26) 
4 no internal fermion loops 
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where G^ 0, is the corresponding 3-point function and S-ij is the quark prop- 
agator. After setting the external gluon momentum equal to zero, the vertex 
can be written as, l8fl: 



A%j(-q,q,0) =gTtj 



7mA 9 (9 2 )+7.(^-^)aJ( 9 2 



(3.27) 



We find the coupling constant that was defined and calculated on the lattice 



in pj, to be: 

gtf) = Avbk (A 9 (-q 2 ) + A T g {-q 2 )) Z^{-q 2 )Z 2 {-q 2 ) , (3.28) 

where 

Sij{q) = -8ij^Z 2 {-q 2 ). (3.29) 
From || one finds, with /x 2 = q 2 in the MS scheme: 

g{q 2 ) = 4^(l + M^]+/i a [^-47C8] + 

, , r 266866067 824999, 349225 , , \ ,„ „„. 
k [ ^7648 28^ C3 -^152- C5] + ---) (3 - 30) 

Putting these coefficients, together with the /3-coefficients (|3.18j) in (|2.5f ) 
(now for N=l/2) and fixing y with the face ( |2.8| ) will give us the two and 
three loop y-expansion results for g(q 2 )- 

The two and three loop MOM scheme results can now be found as 
4vr(/i( n ))^ (n=2,3), with solution of Q . From (gig) and (EOc|) we 
can easily determine the needed /3-coefficients: 

^MOMqg = ™_ 1034C3 ; 
48 

.MOMqg 32456317 4134361 3841475 . 

^ 3 = ~^92" T2 - ^"^SS - Cs • (3 ' 31) 

The A-parameter is now given by: 

2Aln ^M2M3S = ^. (3.32) 

The results are completely similar to the results for the triple gluon vertex. 
Instead of performing a separate fit we simply take the value (237 Mev) 
for Aj^<j obtained from the triple gluon vertex, to compare with the lattice 
result. From fig. |6| we can again observe a turnover for the y-expansion 
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— ■ 2 loop 
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Figure 6: The lattice results |T3 for the coupling from the quark gluon 
vertex with the two and three loop results from the y-expansion and the 
MOM scheme, with A I7? =237 Mev. 



results and the lattice results around 1 Gev, where the MOM scheme results 
diverge. We finally note that the lattice results were obtained for a (small) 
non-zero quark mass, while our results are for a massless quark, so we should 
not be too enthusiastic about the small amplitude difference in the IR-region. 

3.3 The gluon propagator 

The y-expansion will now be applied to the calculation of the scale and 
scheme invariant gluon propagator D^^—q 2 ) defined by: 

d;U~Q 2 ) = f(h)D* (-<?), (3.33) 

with 

» 2 -^2 D U-4 2 ) = (^ + 73^+73^ + . . .)Df v (-q 2 ) (3.34) 

and 

M 2 A/W = ~(73 h + I3,h 2 + 73 2 h 3 + . • .)f(h). (3.35) 
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The general solution of ( 3.35| ) is: 



!{) L i + l/3 ft r + l2/3 2ft + 2/3 ( /3o j 2ft 

(3.36) 



73! _ 73x730/31 73q^1 \,2 

+ 2/?2 /?3 + 2$ r + " 

with A a constant that determines the overall wave function renormalization. 
One can easily check the scale and scheme independence of D. Prom Q and 
(ET36D we find for Z _1 (-g 2 ) (cf. (jlp ), with [i 2 = q 2 and in the WS 
scheme: 

2 v . i, i3„ , r 25085 n l2r 412485993 9747 Al 

Z (-9 ) = \ h~22 (1 + M + 1 r 3 ] + ...). 

V q ' 2904 J L 1874048 352 ^ J ; 

(3.37) 

(Unfortunately we can only determine Z _1 up to second order since for the 
third order result one needs, besides the known third order coefficient for 
Z^ 1 and the four loop /3-coefficient also the four loop 73-coefficient, which is 
not available at the moment. As a consequence we are not able to perform 
an error estimation.) The 2 loop y-expansion result for Z^ 1 is now obtained 
from ( |2.5D , fl2.8| ) and ( 3.3 7| ) . The 2 loop MOM scheme result is found as 
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A 22 where hf^ is the solution of fl2.2|) with 



MOM7 105708585 107217 , 25085 

« ' 29744 20^ Cs and A MOMz = A ¥s ex P^- 

(3.38) 

The 2-loop results for Z(q 2 ) (Euclidean momentum) are shown together 
with a lattice result from |ll| in fig. 0. We now had two fit two things: 
the scale A-j^g and the relative wave function renormalization A. Again, 
we choose the triple gluon vertex value (237 Mev) for Ajfg. A is simply 
determined by fitting the tail of the 2-loop results on the tail of the lattice 
result (at about 5.5 Gev). The overall agreement of our result with the 
lattice is similar as for the vertices. In the deep IR-region, however, there 
is a discrepancy: in it is argued, by extrapolation to infinite lattice- 
volume, that the zero momentum gluon propagator is finite while we find 
a singular zero momentum propagator. Indeed, from the IR-behavior of y 
(jjgj) and the expansion for Z- 1 (p7|) one easily obtains the IR-behavior 
of D(q): 

D(q) ~ 7, — ~ q 01 =q 102 • (3.39) 
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So our zero momentum result is still singular, although the singularity is 
much weaker then the tree level (1/q 2 ) one. We stress that this specific 
power behavior will not be altered by higher loop corrections. 

3.4 The ghost propagator 

The calculation of the ghost propagator is completely similar as for the gluon 
propagator. Again we define the scale and scheme invariant propagator 

G a \q) = -S ab f g {h)G(q 2 ) = -8 ab ^f^. (3.40) 

From |8[ one now arrives at: 

9 1/ 2^ x l, 9 /-, , r 5271 ! ,2r 615512003 5697 . -| 

Z„(-cr) = A" 1 /i"44(i + hi + h 2 l 1 C3 +•••)• 

9 V H ' 9 1936 J 1 7496192 704 ^ ' 

(3.41) 

For the three loop MOM /3-coefficient and the A-parameter we get: 

.MOMgh 653203 6963 . .122 

^ 2 "776 i6~ C3 an A MOMgh = A MS ex P 2904 • 

(3.42) 
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Figure 8: Lattice result taken from fig.l in |l6| (with a 1 =2 Gev) for the 
ghost propagator with the two loop results from the y-expansion and the 
MOM scheme, A I7 cr=237Mev. 



The two loop results for the Euclidean propagator are plotted together with 
the lattice results from 16] in fig. ||. Again we have set Aj^-=237 Mev and 



X g was determined by fitting the two loop results on the lattice results at 
the highest lattice momentum (~ 5.5 Gev, not shown in the fig.). Notice 
the Landau-pole for the MOM result. The agreement of our result with the 
lattice results is satisfying, apart from the strange single data point at the 
lowest lattice momentum. For the IR-behavior of our result we now find: 

, (3.43) 

which is more singular then the gluon propagator but less singular then the 
tree level result. Although this IR-behavior is consistent with |l(|, we should 



remark that other lattice studies [17, Qq] predict a more singular behavior. 



4 Conclusion 

We have presented an alternative perturbative expansion for QCD with only 
one redundant parameter, fixed by a face, and the unexpected feature of 
IR-finite results. The reason behind this was found to be a universal power 
behavior of the running expansion parameter y. For zero flavors there is a 
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qualitative agreement with the lattice data, comparable with the Schwinger- 
Dyson results [19]. 

We do not expect our result to cover all the physics in the intermediate 
energy region and in the IR region. It is immediately clear for example, 
that the y-expansion for non-zero quark flavor will still respect the chiral 
symmetry, so additional non-perturbative corrections are definitely needed. 
However, the fact that our expansion gives finite results in the whole range of 
energies, seems to make it a better (then ordinary perturbation) framework 
to start from, if one wants to estimate the true non-perturbative corrections. 
In the future we will calculate such typical sum-rule quantities as current- 
current 2-point functions. Other possible applications are the calculation of 
experimental cross-sections, starting from the ordinary perturbation theory 
results. 



Acknowledgements 

We would like to thank the lattice community for providing us with some 
of their data. In particular we would like to thank P.Bowman, O.Pene, 
C.Pittori, J.Micheli and J.Skullerud. We also thank K.Chetyrkin for giv- 
ing his (perturbative) data in several output formats. Finally we thank 
O.Bernus for the helping hand on the graphics. 



References. 

[1] P. M. Stevenson. Phys. Rev., D23:2916, 1981. 

[2] G. Grunberg. Phys. Rev., D29:2315, 1984. 

[3] D. V. Shirkov and I. L. Solovtsov. Phys. Rev. Lett., 79:1209-1212, 1997. 

[4] M. Beneke. Phys. Rept, 317:1-142, 1999. 

[5] A. C. Mattingly and P. M. Stevenson. Phys. Rev., D49:437-450, 1994. 

[6] Mikhail A. Shifman, A. I. Vainshtein, and Valentin I. Zakharov. Nucl. 
Phys., B147:385-447, 1979. 

[7] J. Kubo, S. Sakakibara, and P. M. Stevenson. Phys. Rev., D29:1682, 
1984. 



[8] K. G. Chetyrkin and A. Retey. arXiv, [hep-ph/0007088| , 2000. 
[9] P. Boucaud et al. JEEP, 04:006, 2000. 



17 



[10] P. Boucaud, J. P. Leroy, J. Micheli, O. Pene, and C. Roiesnel. JEEP, 
10:017, 1998. 

[11] T. van Ritbergen, J. A. M. Vermaseren, and S. A. Larin. Phys. Lett., 
B400:379-384, 1997. 

[12] William Celmaster and Richard J. Gonsalves. Phys. Rev., D20:1420, 
1979. 

[13] Jonivar Skullerud. Nucl. Phys. Proc. Suppl., 63:242-244, 1998. 

[14] Jonivar Skullerud, Ayse Kizilersu, and Anthony G. Williams. Nucl. 
Phys. Proc. Suppl, 106:841-843, 2002. 

[15] Frederic D. R. Bonnet, Patrick O. Bowman, Derek B. Leinweber, An- 
thony G. Williams, and James M. Zanotti. Phys. Rev., D64:034501, 
2001. 

[16] H. Suman and K. Schilling. Phys. Lett, B373:314-318, 1996. 

[17] Hideo Nakajima, Sadataka Furui, and Azusa Yamaguchi. Nucl. Phys. 
Proc. Suppl, 94:558-561, 2001. 

[18] Attilio Cucchieri, Tereza Mendes, and Daniel Zwanziger. Nucl Phys. 
Proc. Suppl, 106:697-699, 2002. 

[19] Reinhard Alkofer and Lorenz von Smekal. Phys. Rept., 353:281, 2001. 



18 



